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Abstract : The flavor asymmetry of the sea quark distribution as well as the unexpectedly 
small quark spin fraction of the nucleon are two outstanding discoveries recently made in 
the physics of deep-inelastic structure functions. We evaluate here the corresponding quark 
distribution functions within the framework of the chiral quark soliton model, which is an 
effective quark model of baryons maximally incorporating the most important feature of low 
energy QCD, i.e. the chiral symmetry and its spontaneous breakdown. It is shown that 
the model can explain qualitative features of the above-mentioned nucleon structure functions 
within a single framework, thereby disclosing the importance of chiral symmetry in the physics 
of high energy deep-inelastic scatterings. 
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1 Introduction 



Though it is certainly true that high energy deep-inelastic scatterings provide us with the 
best testing ground for perturbative QCD, the physics behind it is highly nonperturbative. In 
fact, the perturbative QCD can describe only the evolution of structure functions from one 
energy scale to the other, so that for determining structure functions themselves one must 
retreat to a semi-theoretical fitting procedure of empirical cross sections at a certain scale [1] . 
In recent years, two epoch-making findings have been made concerning the nucleon structure 
functions. The first concerns the EMC measurement of the spin structure of the proton through 
deep-inelastic scatterings of polarized muons on polarized proton, which sometimes has been 
advertised as the so-called "nucleon spin crisis" [2]. Another interesting observation made by 
the NMC group is the flavor asymmetry of the sea quark distributions in the nucleon [3]. As 
for this second observation, a widely-accepted explanation is based on the pion cloud effects, 
which may in turn be interpreted as a manifestation of chiral symmetry of QCD in high energy 
observables [4-7]. The first observation, i.e. the nucleon spin problem is a little more mysterious 
than the second. The chiral symmetry may play an important role also here [8], or it may be 
explained by some other nonperturbative effects like gluon polarization dictated by the C/a(1) 
anomaly of QCD [9,10]. 

At any rate, these two remarkable findings have reminded us of a quite plain fact : we 
absolutely need some theoretical device with which we can evaluate structure functions in a 
nonperturbative way. Unfortunately, nonperturbative treatment of QCD is quite involved. 
Then, lattice gauge theory is widely believed to be the most promising tool for it [11]. In the 
long run, it may be true. Confining to the present status, however, it suffers from quite a 
few obstacles mainly arising from the limitation of the computer ability. Among others, what 
seems most serious to us is the use of the so-called "quenched approximation", the validity 
of which is strongly suspected when applied to the physics of light quark sector that we are 
just interested in. At least, it seems very unhkely that a rehable estimation of sea quark 
distributions is feasible under this approximation. Here is a place where effective theories of 
QCD can play unique and important roles. Naturally, there exist quite a few effective models 
of baryons, and some of them were already applied for evaluating nucleon structure functions 
[12-14]. Here, we emphasize one absolute superiority of the chiral quark soliton model (CQSM) 
over the other, i.e. its ability to solve the nucleon bound state problem nonperturbatively with 
full inclusion of the deformed Dirac sea quarks in addition to three valence quarks [15-17]. 
Undoubtedly, this must be a unique advantage not shared by other models of baryons when 
discussing the physics of the Gottfried sum [6]. 

Several group have already attempted to calculate nucleon structure functions within the 
CQSM or the NJL soliton model. For instance, Weigel et al. evaluated the isovector un- 
polarized distribution functions as well as some other distribution functions [18]. However, 
their calculation is incomplete in respect that the effect of afore-mentioned Dirac sea quarks 
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are totally neglected. On the other hand, both of Diakonov et al. [19,20] and Tanikawa and 
Saito [21] carried out more consistent calculation including the deformed Dirac sea quarks, 
but by confining themselves to the isosinglet unpolarized and isovector polarized distribution 
functions, which have values at the leading order of 1/A^c expansion (or the expansion in the 
collective rotational velocity Q of the hedgehog soliton). Unfortunately, the isovector unpo- 
larized and isoscalar polarized distribution functions both vanish at this order. Then, for 
discussing interesting physics like the Gottfried sum or the quark spin content of the nucleon, 
one must go ahead to the next-to-leading order contribution in Q, which turns out to be quite 
an involved task. The present paper is the first report of such an investigation. 

The plan of the paper is as follows. In sect. 2, we shall derive general formulas for evaluating 
quark distribution functions of the nucleon on the basis of the path-integral formulation of the 
chiral quark soliton model (CQSM). The results of the numerical calculation are given and 
analyzed in sect. 3. We then summarize our results in sect. 4. 

2 Theory of quark distribution functions 

We start with the light-cone expression for the quark distribution function [12-14] : 

q{x) = - / dz- e-^^''^^- {N{P)\i^{z)Tm\N{P))U^=,,.^=o, (1) 
47r J-oo 

with the definition of the standard light-cone coordinates = {z^ ± z^)/\/2. The nucleon 
state is normalized here as 

{N{P') I N{P)) = 2 P° 5^{P' - P) , (2) 

with P the nucleon 4-momentum, while x = —q^/{2P ■ q) is the Bjorken variable with q the 
4-momentum transfer to the nucleon. We are to take V = 7^X3 with 7"*" = (7° + 7^)/^/2 for 
the isovector unpolarized distribution function u{x) — d{x), while T = 7"^ 75 for the isosinglet 
polarized one Au{x) + Ad{x). We recall the fact that, extending the definition of distribution 
function q{x) to the interval —l<x<l, the relevant antiquark distributions are given as [19] 

u{x) — d{x) — — [u{—x) — d{—x)] , (3) 
Au(x) + Ad{x) = Au(-x) + Ad(-x) , (4) 

with < a; < 1. Taking the nucleon at rest = (Mjv, 0), we have P+ = Mjv/\/2 and then 

q(x) = ^J^Jz'e-'^''-^° {N(P)\^l;\z)^'rm\N(P))Us=-.o,.^^^ (5) 

The basis of our calculation is the following path integral representation of a matrix element 
of an arbitrary nonlocal quark bilinear operator between nucleon states with definite momenta 
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{N{P)\il;\z)OiP{0)\N{P)) = ^ / d^x (fy e-'^-^ e'^y J VU 

X jvil^Vil^^ Jn{^,x) ^\z)0^{Q) Jl{-^,y) exp[z j d^'x^ii^ - MW^)^], (6) 
where 

C = ^^{i ^ - MU^'{x))ilj , (7) 
with W''{x) = exp[i75T ■ 7t{x)/ f^] is the basic lagrangian of the CQSM, and 

is a composite operator carrying the quantum numbers JJ^jTT^ (spin, isospin) of the nu- 
cleon, where is the color index, while ^Yj^ttz^ ^ symmetric matrix in spin-flavor in- 
dices fi. By starting with a stationary pion field configuration of hedgehog shape Uq^{x) — 
exp[i75T ■ rF{r)], the path integral over the pion fields U can be done in a saddle point 
approximation. Next, we consider two important fluctuations around the static configuration, 
i.e. the translational and rotational zero-modes. To treat the translational zero-modes, we 
use an approximate momentum projection procedure of the nucleon state, which amounts to 
integrating over all shift R of the soliton center-of-mass coordinates [24,19] : 

{N{P')\ijj^{z)Oi^{0)\N{P)) — > J d^R {N{P')\ij\zo,z-R)Oij{0,-R)\N{P)) . (9) 

The rotational zero-modes can be treated by introducing a rotating meson field of the form : 

U'''{x,t) = A{t) U;^"{x) A^t) , (10) 

where A{t) is a time-dependent SU (2) matrix in the isospin space. Note first the identity 

^ - MA{t)U^'{x)A\t))ilj = ilj]^{idt- H -n)iljA (H) 

with 

tPa = AUt)jJj , H = + M(3U^'(x) , n = -iAUt)A(t) . (12) 

Here if is a static Dirac Hamiltonian with the background pion fields Uq^{x), playing the role 
of a mean field for quarks, while = ^^aTa is the SU(2)-valued angular velocity matrix later 
to be quantized as Jlo — > Ja/I with I the moment of inertia of the soliton and Ja the angular 
momentum operator. We then introduce a change of quark field variables ip — > ipA, which 
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amounts to getting on a body-fixed rotating frame. Denoting ipA anew as ip for notational 
simplicity, the nucleon matrix element (6) can then be written as 



Z 

X 



{N{P)\ip\z)O^{0)\N{P)) 
J VA Vijj Vi/j^ exp[i J d'^x ^\idt - H - i)' 



Nc i~r i~r Nc rp rri 

X n [A{-)^fX^,x)]ij\zo,z-R)dm-R) n Ki-^,y) ^'i-^)] , as) 

with the definition O = A'''{zo)OA{0). Now performing the path integral over the quark fields, 
we obtain 

{N{P)\ilj\z)Oilj{0)\N{P)) 
^ lfWf{^}^Ar, I d^x d'y e-'P"^ e'Py J d'R J VA 

^ n [/.(o'^l • exp[Splog(za,-//-Q)] , (14) 



j=2 



with r^-^^ = r^-^> )] ' etc. Here Tr is to be taken over spin-flavor indices. Assuming a 
slow rotation of the hedgehog soliton, we can make use of an expansion in VL. For an effective 
action, this gives 

Splog(i9t-//-Q) = Splog(i9t-//) + j nldt. (15) 

The second term here is essentially the action of a rigid rotor, which plays the role of the 
evolution operator in the space of collective coordinates. Using the expansion of the single 
quark propagator as 

f T i i 

- J dz'o dh' /i ( ^ I _ I 4^ Z')a ■ ^a^iz'o) ■ f}{z'o, ^' \ _ jj \ ^0, Z - R)^ 

+ ... , (16) 

we can separate the zeroth order {-^ A^") and the first order (~ 1/A"c) corrections in Q for 
nucleon observables. Here we are interested in the linear order term in Q, because the isovector 
unpolarized and isoscalar polarized distribution functions both vanish at the lowest order. 
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Since necessary manipulations were explained in the previous papers [22,23], we shall skip the 
detail except what seems absolutely necessary. The only but important difference with the 
previous case is that we are here handling a nucleon matrix element of a quark bilinear operator 
which is nonlocal both in space and time coordinates. This is only natural, because we are 
investigating here a quark-quark correlation function with a light-cone separation. First by 
neglecting the time-ordering of two operators Q and O, the O(Jl^) contribution to (14) would 
become 

{N{P)\ij\z)O^{0)\N{P)f' 
= 1 f^^^ f^^>' iVe J d^x d^'y e"*^-^ e'^y j d^R j dh' dz'o J^a/3(4) d^s{zo,0) 



Nc 

X 

J=2 



5 f^^-^l'^'ia;^' • exp[Splog(za,-i/) jnldt], (17) 



As was extensively argued in [22] and [23], one must be careful about the time order of 
two operators VL and O, which do not generally commute after collective quantization of the 
rotational zero-energy modes. In the present case, the proper account of this time-odering can 
be achieved by the following replacement in the equation above : 

^apiZo) d^5{zQ, 0) > [ e{zQ, Zq, 0) + 9{Zq, 0, Zq) ] O^s 

+ [ 0izo, 0, z'o) + 0{O, Zo, 4) ] O^s ^ap , (18) 

where 6 (a, b, c) is a short-hand notation of a step function which is 1 when a > b > c and oth- 
erwise. Here we have intentionally dropped terms containing a factor 6{zo, z'q, 0) or ^(0, Zq, Zq). 
(Note that, if z'q is allowed to lie between zq and 0, the above idea of the time order of the 
two operators Q and O would lose its definite meaning.) It amounts to excluding diagrams in 
which the Coriolis force operates in the time interval between zq and 0. This is motivated 
by the physical picture that a deep inelastic scattering process is a short distance phenomenon 
and its typical time scale is much shorter than that of the collective rotational motion which 
we assume is much slower than the velocity of the intrinsic quark motion. Then, one should 
rather treat Zq and as nearly degenerate when performing the expansion in Q, which dictates 
the neglect of diagrams containing 9{zo,Zq,0) or 9{0,z'q,Zo). The exclusion of these special 
time-order diagrams has an important physical consequence as we shall come back later. An- 
other important remark is that we do not give any constraint on the time order of zo and 
themselves : zq can be earlier or later than 0. It can be verified that inclusion of both 
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possibilities is essential for maintaining the charge conjugation symmetry, which denotes that 
the quark (antiquark) distributions in the nucleon must be the same as the antiquark (quark) 
distributions in the antinucleon. The last comment we want to make here is on the treatment 
of the rotated operator O — A\zq) 70 O A(0). If the operator contains an isospin factor (or it 
is an isovector operator), we are to use the identity 

A^{zo)TaA{Q) = ]^Tx{A\zo)TaA{Q)n)n. (19) 

A new feature here is that the factor Tr(A^(2;o) A(0) r^) contains two time arguments zq 
and 0. Still, we can use the usual procedure in which it is replaced by the standard Wigner 
function Dah{A). Intuitively, the justification follows from our dynamical assumption that zq 
and can be regarded as nearly degenerate as compared with the time scale of the collective 
rotational motion. A formally more rigorous justification for this replacement may be obtained 
by remembering that the time evolution of the collective space operator A is given by 

A{zq) = exp(izo ^7^) ^(0) exp(-izo JV^) > (20) 

with J the angular momentum operator working in the collective coordinate space and / the 
moment of inertia of the soliton. Since the moment of inertia / is an 0{Nc) quantity, this gives 

^ IV {A\zo) Ta ^(0) n) = Dab{A) [1 + 0{l/N,)] , (21) 

convincing that the error of the above replacement is a higher order correction in the 1/A^c 
expansion. 

After stating all the delicacies inherent in the structure function problem, we can now 
proceed in the same way as [22] and [23]. Using the spectral representation of the single quark 
Green function 

a(a;,t|— ^|a;',t% = e{t-t') e-''^-^'-'K{x\m){m\x')p 

^Ot - n m>Q 

- e{t' -t) e-^^-(*-*') o,{x \m){m\ x')p , (22) 

m<0 

together with the relation 

{z-R\ = {-R\e'P-^ (23) 

we can perform the integration over R, z', and Zq. The resultant expression is then put into (5) 
to carry out the integration over zq- We then arrive at the formula, which gives the theoretical 
basis for evaluating the first order contributions in Q, to quark distribution functions of the 
nucleon : 

q{x;n') - J dA ¥j%[A] O^'^A] ¥fJ^^[A] . (24) 
Here 



= . (25) 



are wave functions, describing the collective rotational motion of the hedgehog soliton, while 



0^'\A] = + O^liA] , (26) 

where 

^ { !L (0 I O \m) {m\fl\0) 

+ I O z I 0) (0 1 1] I m) } 

+ (O^Q) . (27) 



and 



X { < n I O I m) {m\\l\n) + {O \l) \ , (28) 

with the ^-component of the single-quark momentum operator. In the above formulas, | m) 
and Em denote the cigcnstatcs and associated cigenenergies of the static Dirac Hamiltonian 
H . In particular, | 0) represents the lowest energy eigenstate of iJ, which emerges from the 
positive energy Dirac continuum. This particular state is referred to as the valence quark 
orbital following the convention of previous studies [16-17]. However, the term "valence quark" 
here should not be confused with the corresponding term in the language of the quark parton 
model. The two contributions 0^^\ and O^^^ can actually be put together into the form : 



0« = M^.N, 



E — E 

m>0,n<0 

x{(„|6^'"^"-^--'''>t^'^^"-^--''''|m)(,»|a|») + (6»0)}. (29) 

which contains transitions from occupied single quark levels to unoccupied ones in consistent 
with the Pauli principle [23]. It is important to recognize that this nice feature would not 
follow, if one would include the afore-mentioned particular time-order diagrams containing 
d{zi^, z'q, 0) and ^^(0, z'q, zq). In fact, if we had included these contributions by forgetting about 
the time order of Q and O, we would have obtained 



2 ^ E — E 



, / I /s sign(£^m) S{xMn -Em- p^) - sign(K) S{xMn - E^ - p^) . , , , „ , , 

X (n I O ^ — — ^ — - — \m) {m\il\n) , (30) 



instead of (28). This essentially coincides with the equation which is obtained from (4.12) 
in ref.[19]. Clearly, the formula (30) contains unpleasant Pauli- violating contributions due to 
transitions between occupied levels themselves and unoccupied ones. It is very important to 
recognize that this Pauli- violating contributions to a quark distribution function has nothing 
to do with the regularization problem. (This is clear, since we have not yet introduced any 
regularization at this stage.) This kind of Pauh-principle violation has never been observed 
when only the nucleon matrix elements of local quark bilinear operators have been dealt with. 
In fact, one can confirm it very easily, if one carries out the x integration of (30) from x — —Ito 
X = 1, which should give the first moment of the relevant distribution function being expressed 
as a nucleon matrix element of a local quark bilinear operator. One certainly convinces that 
terms with x > drop completely. 

We also call attention to the fact that the delta-function representing the on-shell energy 
conservation appears in the symmetrized form with respect to the energies of the single-quark 
levels I m) and | n). The appearance of these two delta- functions can be traced back to two 
possible time order of Zq and in 'ip{z)Tip{0), and it plays a favorable role to maintain the 
charge-conjugation symmetry as already mentioned. We point out that a naive cranking 
procedure without taking care of the non-local (in time space) nature of the quark bilinear 
operator would violate this property [18]. 

Using the general formulae given in (24) ~ (28), we can readily obtain the following ex- 
pression for the isovector unpolarized distribution function for the proton : 



(In the actual numerical calculation, we use the expression in which the contribution of the 
valence quark level is separated from that of the Dirac continuum according to the general 
formulas (26) ~ (28), since only the vacuum polarization part is to be regularized in our 
regularization scheme.) 

By integrating u{x) — d{x) over x between —1 and 1, the isospin sum rule follows; 




a=l ^ m>0,n<0 "^"^ 

2 



X {n\Ta{l + 7° 7') 



m) (m I Ta I n) . 



(31) 




{ [u{x) — d{x)] — [u{x) — d{x) ] } dx 



(32) 



J-i Jo 



provided that the moment of inertia is given by 




(33) 



On the other hand, the Gottfried sum is given as 




(34) 
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The formula for the isosinglet polarized distribution function can similarly be written down 

as 



IN 1 
Au{x) + Ad{x) = -M^.-.^ J: 

^ m>0,n<0 ^■rn 

X (n 1(1 + 7''7'') 75 I m) (m I Ts I n) . (35) 

Integrating over x, we get 

j [Au{x) + Ad(x) ] dx ^ J [Au{x) + Ad{x) + Au{x) + Ad{x) ] dx 

= -y--^ F i F (^|7°7^75|m) (m|r3|n) , (36) 

^ m>0,n<0 -^"^ 

which is noting but the isosinglet axial- vector charge ( or the quark spin fraction) of the nucleon 
first investigated in [16] in the context of the CQSM. 



3 Numerical results and discussion 

Now we turn to the discussion on the method of numerical calculation. For evaluating 
isovector unpolarized and isoscalar polarized distribution functions with inclusion of the Dirac- 
sea quark degrees of freedom, we must perform infinite double sums over all the single-quark 
orbitals which are eigenstates of the static hamiltonian H. A numerical technique for carrying 
out such double sums has been established in the case of static nucleon observables [16]. In 
the nucleon structure function problem, however, we encounter a new bothering feature. The 
relevant matrix elements contain delta functions depending on the single-quark momentum 
operator p^. Kahana and Ripka's plane wave basis [25] is useful also for treating this feature. 
(We recall that it is a set of eigenstates of the free hamiltonian Hq = + (3M discretized 
by imposing an appropriate boundary condition for the radial wave functions at the radius 
D chosen to be sufficiently larger than the soliton size.) For evaluating necessary matrix 
elements, we find it convenient to expand the eigenstates | m) and \ n) oi H in terms of this 
discretized plane wave basis in the momentum representation. We can then make use of the 
fact that is diagonal in this basis. This however does not get rid of all the trouble. Since the 
delta function in question now depends on the Bjorken variable x, Em (or £■„), — \p\i cosOp, 
where \p\i are the momenta of the discretized basis states, the resultant distribution would be a 
discontinuous function of x [20]. To remedy it, Diakonov et al. proposed a smearing procedure 
in which a calculated distribution function is convoluted with a Gaussian distribution with 
an appropriate width [20]. This smearing procedure works well as a whole, but it turns out 
to have some unpleasant features. First, it requires us a special care near a; ~ 0, where 
a distribution function (extended to the range — 1 < x < 1) or its derivative may have a 
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discontinuity. Secondly, the smearing causes unpleasant flattening when the bare distribution 
has a sharp peak as is the case for the valence quark contribution. Then, we decided to use 
here a conceptually simpler least square fitting procedure. To this end, we first calculate 
distribution functions at many points of x with a small interval, say Ax = 0.01. This gives a 
set of data, which are rapidly fluctuating with a small amplitude. A common variance of the 
order of this fluctuation amplitude is then given to all the data points, so that the least-square 
fit can be carried out. In the fitting of the distribution functions, wc find it convenient to use 
a set of locally peaked Gaussian functions, the centers of which are distributed homogeneously 
in the x space. 

Actually, the expression for the isovcctor unpolarizcd distribution as well as that for the 
moment of inertia contain divergences, so that they need some regularization. In the case of 
static nucleon observablcs, there is a phcnomcnologically established regularization method, 
based on Schwinger's proper-time scheme [16,17]. Unfortunately, how to generalize this reg- 
ularization scheme in the evaluation of nucleon structure functions has not been known yet. 
In the recent investigations of the leading order quark distribution functions [19,20], Diakonov 
et al. advocated to use the Pauli-Villars regularization. In this scheme, one subtracts from 
divergent sums a multiple of the corresponding sums over eigenstates of the hamiltonian in 
which the quark mass M is replaced by the regulator mass Mpv, uniquely determined from 
the condition, 

- 4^2 ' ^"^^^ 

which gives Mpy ~ 571 MeV for M = 400 MeV. In the present study, we shall also try this 
regularization scheme but somewhat in a different way. That is, although they regularized 
not only the Dirac continuum contribution but also the discrete (valence) level one, it is not 
an absolute demand of this regularization scheme. Since the valence level contribution is 
convergent itself, one can also define and use another regularization scheme, in which only the 
Dirac sea part is regularized by the Pauli-Villars subtraction. In fact, a closer correspondence 
with the standardly-used regularization scheme like the proper-time one is held by adopting 
this latter scheme of regularization. Either way, since this method of regularization with use 
of subtraction method is quite dissimilar to more standard energy-cutoff scheme, one may 
suspect that the answers would depend strongly on this particular choice of regularization. To 
see how much the predicted distributions depend on the choice of the regularization method, 
we have invented another semi-theoretical regularization scheme, which we believe is closer in 
spirit to the standard method of introducing regularization. Let us explain it for u{x) — d{x), 
since the vacuum quark part of /\u{x) + Ad{x) is convergent without regularization. We start 
with the nonregularized form for the vacuum polarization part of u{x) — d{x) 

a=l ^ m>0,n<0 
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X < n I Ta ( 1 + 7 7 ) ^ — — ^ \ m> < m\T3\n> . (38) 

One of the simplest choice for introducing a regularization, which does not violate the charge- 
conjugation symmetry, would be given by the replacement, 

b{xM^ - K - P^) + K^M^ -Em- p^) 



2 

g{Er,; E^ax) S{xMn -E^- p^) + g{Em; Emax)S{xMN - E^ - p^) 



(39) 



2 

where g{En] Emax) is an appropriate energy cutoff which may, for instance, be taken as 

g{Er,]Ernax) = Crfc (| ^ | / ^max ) , (40) 

or 

g{En,E„,ax) = eM-El/Elax)- (41) 

For definiteness, we will take the second choice in the following explanation. Integrating over 
X from —1 to 1, this leads to the Dirac sea contribution to the moment of inertia of the form : 



ax 



Iv.p. ^ ^ F _ F o {n\T3\m) {mlnln), (42) 

^ m>0,n<0 ^■"T' ^ 

where use has been made of the symmetry of the matrix element of t^. One notices that 
this expression is pretty different from that of the standardly-used proper-time regularization, 
which gives 

Iv.p. = 13 f{Em, En] A) < n I Ta I m >< m I Ta I n >, (43) 

m,n 

with 

P 1 , sign(E^) erfc(|E^|/A) - sign(EO erfc(|K|/A) 

J [Em: En -.A) = T "I ^ ^ 

-T^' El-El 
though both are reduced to the following same expression in the infinite cutoff hmit : 

Nc 1 

Iv.p. = ^ XI _ < n\T2,\m >< m\T3\n > . (45) 

^ m>0,n<0 -^"i 

A noticeable difference is that the formulas based on the proper-time regularization contains 
transitions between occupied levels themselves and unoccupied ones, which clearly violates 
Pauli exclusion principle. However, this physically unacceptable Pauli-principle violation of 
the proper-time regularization is a mere appearance. Due to the favorable function of the 
second term of the cutoff function f {Em, En, A), it turns out that the contribution of the 
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Pauli-violating processes with x En > turns out to be greatly suppressed. On the other 
hand, the regularized expression (42) is free from the Pauli-principle violation by construction. 
However, since the behavior of the cutoff functions in (42) and (43) are pretty different, there 
is no reason to expect that the cutoff parameter E^ax in (41) is close to A in (44). That means 
that we need some criterion to determine the parameter E^ax in (40) or (41). In view of the 
fundamental importance of the moment of inertia in the collective quantization approach, a 
natural procedure would be to fix the cutoff parameter so as to reproduce this key quantity. 
One may then fix Ej^ax so that the new regularization scheme gives the same moment of inertia 
as the standard proper-time scheme. Here, for the sake of convenience, we would rather fix it 
so that the new regularization scheme gives the same moment of inertia as the Pauli-Villars 
regularization scheme. The reason is as follows. In the previous studies of static nucleon 
obscrvablcs, all the observables can be evaluated within the single regularization scheme, which 
also enables us to get a self-consistent soliton solution [16,17]. In the present semi-theoretical 
treatment of the regularization problem, we must give up this complete self-consistency. Under 
this circumstance, the best starting choice would be to use a phenomenologically tested soliton 
profile, which has been obtained in the proper-time regularization scheme with use of the 
dynamical quark mass M — 400 MeV and the physical pion mass = 138 MeV. Using this 
soliton profile, we obtain Ifotai — 1.16 fm in the proper-time regularization scheme, while the 
Pauh-Villars regularization scheme with the same profile function gives Itotai — 1.14 fm. This 
means that the above two criterions for fixing E^ax can be thought of as practically the same. 
Now by setting as 

E„,ax = 728 MeV for (?(E„;E„„,) = exp{-El/ E^^J , 

En^ax = 1080 MeV for g{En;E„,ax) = eric{\En\/E^ax) , (46) 

the new regularization schemes give 

/^.p. (Gaussian-function) ~ 0.254 fm , 

(error- function) ~ 0.254 fm (47) 

which approximately reproduce the number 

7„.p. (Pauli-Villars) ~ 0.255 fm, (48) 

obtained with the Pauli-Villars regularization scheme. Thus, at least once an appropriate 
soliton profile is given, cither of these three regularization schemes thought expected to be a 
reasonable candidate for regularizing divergences in our problem. In the following, we shall 
calculate the vacuum polarization contributions to u{x) — d{x) and u{x) — d{x) by using these 
three different schemes of regularization, to see the regularization dependence of distribution 
functions,. 
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Now we are in a position to show the results of our numerical calculation. Shown in Fig.l 
are the contributions of the (discrete) valence level to the distribution functions u{x) — d{x) 
and u{x) — d{x) in the proton. One clearly sees that u{x) — d{x) is negative in the range 
of X where it has dominant support. This means that the excess of the d sea over the u in 
the proton comes out very naturally even at the level of valence quark approximation. This 
feature, which was also observed in ref.[18], is an interesting consequence of the CQSM, which 
assumes the symmetry breaking mean-field of hedgehog form. One may be interested in what 
conclusion one would obtain if one evaluates the net flavor asymmetry of sea quarks under this 
valence quark approximation. One gets 

/' [uix) - dix)l^^ dx ^ -0.017. (49) 

t/ 

(Here the upper limit of x integration is actually extended to oo, since the theoretical distri- 
bution functions have non-zero support beyond x — 1. See the discussion before (52).) The 
above number may be compared with the empirical estimate 

[\u(x) - d(x)] dx ~ -0.148, (50) 
Jo 

which is obtained from the NMC measurement Sq = 0.235 ± 0.026 [26]. One sees that the 
valence quark contribution alone is not enough to generate the required magnitude of flavor 
asymmetry of the sea quark distributions. As a matter of course, a comparison with the 
experimental data is premature at this stage of calculation. First of all, we have not yet paid 
any attention to the renormalization scale dependence of the distribution functions. (We shall 
later come back to this point.) Secondly, the distribution functions of Fig.l obtained within 
the valence quark approximation do not saturate the isospin (or Adler) sum rule. In fact, we 
obtain 

/' { [u{x) - d{x) ] - [u{x) - d{x) ] dx = 0.775 , (51) 

J 

denoting that the remaining 23% of the nucleon isospin is carried by the Dirac sea quarks. 
The isospin sum rule can be made to hold within the valence quark approximation if the total 
moment of inertia / in (31) is replaced by its valence quark contribution I^ai, as was done in 
[18]. However, this would also change all the O(Q^) distribution functions by the factor of 
1 1 lyai — 1.3, which is by no means justiflable. 

Anyhow, it is clear that a consistent calculation should also include the Dirac sea contri- 
bution. Since this part is generally dependent on the selected regularization scheme, we first 
investigate in Fig. 2 the regularization dependence of the vacuum polarization contribution 
to u{x) — d{x) and u{x) — d{x). Here, the solid curves are the results of the Pauli-Villars 
regularization scheme, while the dashed and dash-dotted curves arc respectively those of the 
energy-cutoff scheme with the Gaussian and error-function type regularization functions. One 
sees that the regularization scheme with use of the energy cutoff has a tendency to suppress 
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distribution functions at small x as compared with the Pauli-Villars scheme. However, the 
differences between the three regularization scheme are not so large. One can rather say that 
the regularization scheme dependence of the distribution functions is pretty small once the cut- 
off parameters are determined so as to reproduce a reasonable value of the moment of inertia. 
Note also that the regularization scheme dependence becomes further insignificant if we see the 
sum of the valence and the vacuum polarization contributions to the structure functions, since 
the contribution of the valence level is dominant anyway. Shown in Fig. 3 arc the sums of the 
valence and the vacuum polarization contributions to u{x) — d{x) and u{x) — d{x). Here only 
the distribution functions obtained with the Gaussian-type regularization function is shown 
for the reason explained above. We have numerically confirmed that the isospin sum rule (32) 
is satisfied within the precision of 0.2 %. One sees that, as compared with the distribution 
functions obtained with the valence quark approximation (Fig.l), a further enhancement is 
observed for the excess of the d sea over the u sea inside the proton. 

So far we have postponed the discussion on the energy scale of the model calculation. 
Roughly speaking, the quark distribution functions evaluated here corresponds to the energy 
scale of the order of the Pauli-Villars mass Mp.y. — 600 MeV or the cutoff energy Ej^ax — 
700 MeV contained in the Gaussian-form regularization function. The Q^-evolution must be 
taken into account in some way to be compared with the observed nucleon structure functions 
at large Q'^. Fortunately, after elaborate analyses of high data by taking account of the 
perturbative evolution, Gliick, Reya and Vogt gave a simple parameterization of quark 
distribution functions at a normalization point very close to the energy scale of the present 
model [29]. We may therefore carry out a preliminary comparison of our predictions with 
their low scale parameterization of quark distribution functions. Fig. 4 shows this comparison. 
Here the solid curves stand for the distribution functions x [u{x) — d{x) ] and x [u{x) — d{x) ] 
obtained from the ones given in Fig. 3, whereas the boxes represent the NLO parameterization 
of ref. [29]. One can say that the qualitative feature of the NLO parameterization is nicely 
reproduced. The main discrepancy between the theory and the phenomenological fit is that 
the theoretical distribution functions have nonzero support beyond x — 1. This unphysical tail 
of the theoretical distribution functions comes from an approximate nature of our treatment 
of the soliton center-of-mass motion (as well as the collective rotational motion), which is 
essentially nonrelativistic. A simple procedure to remedy this defect was proposed by Jaffe 
based on the 1-1-1 dimensional bag model [27] and recently reinvestigated by Gamberg et al. 
within the context of the NJL soliton model [28] . (The latter model is essentially equivalent 
to the CQSM.) According to the latter authors, the effect of Lorentz contraction can simply 
be taken into account by first evaluating the distribution functions in the soliton rest frame 
(as we are doing here) and then by using the simple transformation 



as far as the order O contributions to the distribution functions are concerned. We are not sure 
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9(1 -x) 



(52) 



whether their proof can be generahzed for the order fl^ contributions to distribution functions, 
in which we must treat three dimensional collective rotational motion. Nonetheless, it may be 
interesting to see the effects of this transformation, especially on u{x) — d{x). The distribution 
functions obtained after this transformation are shown in Fig. 4 by the dashed curves. The 
crucial influence of this transformation at large x region is obvious from this figure. Since (52) 
is a normalization-preserving transformation, i.e. 

/ fRF{x) dx ^ fiMF{x) dx , (53) 
Jo Jo 

it follows that the peaks of the transformed distribution functions becomes sharper than the 
original ones. However, the behavior of the distribution functions at smaller values of x turns 
out to be rather insensitive to this transformation. In this small x region, there remains a 
qualitative difference between the phenomenological distribution functions and the theoretical 
ones (obtained after the transformation (52)) especially for the antiquark distribution u{x) — 
d{x) . Note however that the enhanced behavior of the GRV parameterization at smaller values 
of X is connected with the Regge-like behavior (~ l/x" with < a < 1) of the distribution 
functions assumed in their fit. It is clear that such hidden dynamics of the distribution functions 
never enters into the model calculations as carried out here. Probably, within the model 
calculation as carried out here, what is more reliable than the detailed x dependence would be 
sum rules obtained after the x integration. Integrating over x, we obtain from Fig. 4, 

5g = [\u(x) - d(x)] dx ~ -0.130, (54) 
Jo 

or equivalently 

Sg ^ - [\[u(x) - d(x)] + [u(x) - d(x)]} dx ~ 0.247. (55) 
3 Jo 

On the other hand, using the NLO parameterization of ref.[29], we find that 

5g ^ - 0.168 , (56) 

or 

Sg ^ 0.221 . (57) 

We can say that the agreement between the theory and the phenomenological fit is pretty 
good in view of the extreme simplicity of the CQSM. Needless to say, the flavor asymmetric 
polarization of the sea quark distributions obtained above would never arise if there is no flavor 
asymmetry in the valence quark numbers of the nucleon. This then denotes that the NMC 
observation is nothing mysterious : it is explained very naturally as a combined effects of two 
ingredients, i.e. the apparently existing flavor asymmetry of the valence quark numbers in the 
nucleon and the spontaneous chiral symmetry breaking of the QCD vacuum [30] . 
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Next we turn to the discussion on the flavor-singlet (isoscalar) polarized distribution func- 
tion. As has been already pointed out, the vacuum polarization contribution to this quantity 
is convergent without regularization. (This feature comes from the fact that the flavor-singlet 
axial-charge is related to the imaginary part of the Euchdean effective meson action.) What is 
more, we flnd that the Dirac sea contributions to Au{x) + Ad{x) as well as Au{x) + Ad{x) are 
consistent with zero within the numerical accuracy of the model calculation. This means that 
the valence quark approximation is a good approximation for this special quantity. However, 
the valence quark approximation here should be distinguished from the corresponding one in 
ref.[18] where the total moment of inertia / had to be replaced by its valence quark part I^ai 
so as to maintain the isospin (or Adler) sum rule. In making this replacement, one would 
necessarily overestimate 0{Q^) distribution functions like the isoscalar polarized one, since 
they are inversely proportional to /. In any case, since the vacuum polarization contribution 
to Au{x) + Ad{x) as well as Au{x) + Ad{x) are almost negligible, we will show below only 
the total contributions. The solid curves in Fig. 5 stand for the isoscalar polarized distribution 
functions Au{x) + Ad{x) and Au{x) + d{x) directly obtained from the formula (35), whereas 
the dashed curves here take account of the effect of Lorentz contraction by making use of the 
transformation (52). Again, the sizable effect of Lorentz boost is self-explanatory especially 
at larger values of x. Since the low energy scale parameterization of the polarized distribution 
functions is not yet available, we shall postpone the detailed comparison with the observed 
structure functions for future studies. A qualitatively interesting feature observed in Fig.5 is 
that the isoscalar combination of u and d seas is slightly polarized (in most range of x) against 
the nucleon spin direction. Although this combined polarization of u and d seas is rather small 
in magnitude, this does not mean that Au{x) and Ad{x) are both small. The leading order 
calculation of Diakonov et al. for the isovector polarized distribution functions tells us that 
Au{x) — Ad{x) is likely to be positive and fairly large [19,20]. This together with the above 
result for Au{x) + Ad{x) indicates that the seas of u and d are sizably polarized in such a 
way that their polarizations cancel each other. Now using the distribution functions given in 
Fig.5, we can calculate the first moment of the fiavor singlet (longitudinally) polarized quark 
distribution functions. We flnd that 

C [ Auix) + Adix) + Auix) + Adix) ] dx = 0.494 . (58) 

On the other hand, the direct calculation of the flavor-singlet axial charge with use of (36) 
gives 

{gf) = 0.492. (59) 

The agreement between the above two numbers can be interpreted as showing the accuracy 
of our numerical procedure for evaluating quark distribution functions. In any case, from the 
above value of the sum rule, one reconflrms that the spin fraction of the nucleon carried by 
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quarks comes out to be quite small (less than half) in the CQSM. As has been repeatedly em- 
phasized by one of the authors [16], the smallness of the quark spin fraction (or the largeness of 
the orbital angular momentum) is inseparably connected with the basic dynamical assumption 
of this unique model, i.e. the identification of a rotating hedgehog with the physical nucleon. 

4 Summary 

In summary, we have evaluated the isovector unpolarized and isoscalar polarized distribution 
functions of the miclcon within the framework of the CQSM with full inclusion of the Dirac 
sea quarks . It has been shown that the characteristic features of the observed distribution 
functions, i.e. the excess of the d sea over the u sea in the proton as well as the very small quark 
spin fraction of the nucleon, are reproduced at least qualitatively within a single theoretical 
framework of the CQSM, which is the simplest effective quark model of QCD taking maximal 
account of chiral symmetry. We then interpret this success as revealing the crucial importance 
of chiral symmetry in the physics of high energy deep-inelastic structure functions. 
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Figure caption 



Fig. 1. The contributions of the valence level to the distribution functions u{x) — 
d{x) and u{x) — d{x). 

Fig. 2. The contributions of the Dirac continuum to the distribution functions 
u{x) — d{x) and u{x) — d{x). The solid curves represent the results of 
the Pauli-Villars regularization scheme, whereas the dashed and dash- 
dotted curves are those of the energy-cutoff scheme with the Gaussian 
and error-function type regularization functions. 

Fig. 3. The sum of the valence and Dirac sea contributions to u{x) — d{x) and 
u{x) — d{x) obtained with the regularization function of Gaussian form. 

Fig. 4. The theoretical distribution functions in comparison with the low energy 
scale parameterization (NLO fitting) of ref. [29]. The solid curves just 
correspond to the distribution functions of Fig. 3 multiplied by x, while 
the dashed curves are obtained by using the transformation (52). 

Fig. 5. The isoscalar polarized distribution functions AM(x)+A(i(a;) and A-u(a;) + 
Ad{x). The solid curves are directly obtained from (35), while the dashed 
curves are obtained with the transformation (52). 
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